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Iron sesquioxide (Fe20s) displays pressure and temperature induced spin and structural transi- 
tions. Our calculations show that, density functional theory (DFT), in the generalized gradient 
approximation (GGA) scheme, is capable of capturing both the transitions. The ambient pressure 
corundum type phase (hematite or a-Fe203), having R% c symmetry, gets distorted by the appli- 
cation of pressure and transforms to a distorted corundum type or Rh20s(H) phase with Pbcn 
symmetry, in agreement with recent experiments. GGA + U calculations show the same trend but 
shift the transition pressures to higher values. Experimentally, the onset of the structural transition 
begins in the vicinity of the spin transition pressure and whether the system undergoes spin tran- 
sition in the corundum type (HP1) or in the Rh20s(II) type (HP2) phase, is still a controversial 
issue. With a relatively simple, but general, octahedral structural parameter, V oc t (the octahedral 
volume around iron ions), we show that in order to acquire a low spin (LS) state from a high spin 
(HS) one, the system does not necessarily need to change the crystal structure. Rather, the spin 
transition is a phenomenon that concerns the cation octahedra and the spin state of the system 
depends mainly on the value of V oc t, which is governed by two distinct equations of state, separated 
by a well defined volume gap, for the HS and LS states respectively. Analysis of the results on the 
basis of octahedral volume allows to sum up and bridge the gap between two experimental results 
and thus provides a better description of the system in the region of interest. 

PACS numbers: 



I. INTRODUCTION 

Spin crossover or electronic transition is a phe- 
nomenon commonly encountered among (transition) 
metal complexes, in particular under octahedral 
coordination 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 . Under the influence of the ligand 
field the metal d orbitals are split into two sets: a tig level 
with a 3-fold degeneracy and a 2-fold degenerate e g level. 
The closer the ligand is to the metal ion, the greater is the 
influence of the ligand-field strength. When this ligand- 
field energy dominates over electron-pairing or exchange 
energy the metal ion adopts the low-spin configuration. 
A high-spin configuration is instead favored when the ex- 
change energy is greater than the ligand-field energy. 

Hematite, a-Fe203, is a rhombohedrally structured 
wide-gap antiferromagnetic insulator 8 at ambient condi- 
tions. It becomes weakly ferromagnetic between Morin 9 
temperature, Tm (=260 K), and Neel temperature, Tn 
(= 955 K), as a result of the canting of the spins of the two 
sublattices 10 ' 11 ' 12 ' 13 . In this structure, the Fe 3+ ions, sit- 
ting in octahedral cages, display a high spin (HS) to low 
spin (LS) transition under external influences such as, 
temperature, pressure, etc. Sequential or concomitant 
pressure and/or temperature induced structural phase 
transitions, in addition to the spin one, increase the com- 
plexity of the system, making its in-depth understanding 
a challenge. 

Significant continuous attention by the scientific com- 
munity, for more than forty years, has brought in many 
intuitive ideas and counter-ideas 14 ' 15 ' 16 ' 17 . The existence 
of residual magnetism after a volume collapse of ~ 10 % 
and concurrent onset of the structural phase transition 
led to the issue of one-type or two-type cationic picture 



for the high pressure phase 18 ' 19 ' 20 ' 21 . To resolve this is- 
sue and correlate/separate the spin and structural tran- 
sitions a number of investigations have been performed 
in more recent years. Of them, X-ray diffraction (XRD), 
Mossbauer spectroscopy (MS) at 300K and electrical re- 
sistivity measurements by Pasternak et. al. 22 assigned 
a nonmagnetic distorted corundum or RI12 03(H) type 
structure as the high pressure phase. The intermediate 
region between the insulating and metallic regions in the 
electrical resistivity data was qualitatively explained by 
the coexistence of an insulating magnetic phase with a 
metallic nonmagnetic one. 

Another XRD and X-ray emission spectroscopy (XES) 
experiment by Badro et. al 23 could separate the elec- 
tronic transition from the crystallographic one. From 
the XES data these authors concluded that a high-spin 
state could be stabilized in the high-pressure structural 
phase at low temperatures while at high temperatures 
a low-spin state was stable. This high temperature LS 
state was found to be still weakly paramagnetic as evi- 
denced by the presence of a satellite (K/?') in the XES 
spectra 23 , in contrast with the non-magnetic state sug- 
gested by the authors of Ref. 22. The structural tran- 
sition shows a sluggish behavior, as reported in room 
temperature Raman spectroscopy studies by Shim and 
Duffy 24 and angle-resolved XRD studies by Rozenberg 
et. al. 25 . Moreover, high temperature heating might have 
a role to play in controlling the transition pressure as 
shown in the XRD experiments by Ono et. al 26,27 . 

Hence, in spite of significant experimental efforts de- 
voted to this system, some still controversial and unre- 
solved issues remain: (i) how does the spin state of the 
system evolve with pressure ? (ii) what is the correlation 
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between the spin and the structural transition, if there 
is any? In this study we will show that spin transition 
is a general phenomenon, depending on a simple struc- 
tural parameter, the octahedral volume around Fe ions, 
and thus largely uncorrelated to the structural transition 
occurring in the same pressure range. The apparent con- 
tradiction between experimental results reported in Ref. 
22 and 23 will also be resolved. 



II. METHODS AND RESULTS 

Extensive first-principles density functional calcula- 
tions of the spin and structural transitions in Fe203 
have been performed. Variable-cell-shape molecular 
dynamics 28 has been used for the full structural optimiza- 
tion at arbitrary pressures. Eight-electron (3d 7 4s 1 ) and 
six-electron (2s 2 2p 4 ) ultrasoft pseudopotentials 29 wprp 
used to describe Fe and O atoms 30 . The electron- 
correlation energy was treated within the PBE 1 
alized gradient approximation (GGA). GGA+U 
tions were performed to supplement and comp 
the results obtained within GGA. In these Cc 
U values have been calculated self-consistentl 
the methodology described in Ref. [32]. All cal< 
presented in this work have been performed 1 
Quantum-ESPRESSO 33 package. 

To understand the evolution of the spin and si 
transition in the system we have considered tf. 
sible structural phases, starting from ambient 
pressure. Fig. 1(a) refer to the rhombohedral 
tured corundum- type Fe2 03, where the four I 
are arranged linearly along the rhombohedral a> 
unit cell. The two possible candidate structur 
evated pressure are Rli203(II)-type and orthoi 
perovskite (opv) type Fe203, as shown respect 
Fig. 1(b) and Fig. 1(c) and contain four formi 
per unit cell. For each of the above mentioned st 
there can be a few anti-ferromagnetic conflguri 
addition to the ferromagnetic one, depending on 
orientations. 




FIG. 1: (Color online) Structures for the (a): corundum- type, 
(b): Rh20s(II)-type and (c): optotype Fe2C>3. 



To start with we present the density of states (DOS) 
for each of the abovesaid three structures for their most 
stable ambient-pressure magnetic state. As can be seen 
from the lower panel of Fig. 2, DFT at the GGA level cor- 
rectly describes the insulating nature of hematite. How- 
ever, there is, as typical of LDA/GGA calculations, an 
underestimation of the band gap. Rli203(II)-type Fe203 
also shows the same behavior, but with a much dimin- 
ished magnitude of the band gap, as shown in the middle 
panel of Fig. 2. In both cases the valence band top edge 
is equally populated by Fe-3<i and 0-2p states. However, 
DOS for the opv-type Fe203, shown in the upper panel 
of Fig. 2, displays a very poor metallic tendency with 
Fe-3d dominated valence edge. This could be due to the 
underestimation of correlation effects by GGA. However, 
this issue is not relevant here since the latter structure is 
never the most stable one in the pressure range of inter- 
est, as will be shown below. 
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FIG. 2: (Color online) Projected density of states calculated 
at GPa within GGA. Lower panel: corundum- type. Mid- 
dle panel: Rh20s(H)-type and upper panel: opv-type Fe203- 
Note that in the upper panel the dotted lines (Fel -3d) cor- 
respond to Fe atoms with 8-fold coordination. 

The calculated relative enthalpy curves of Fe203 for 
corundum- type, Rli2 03(II)-type and opv-type structures 
in different magnetic configurations are reported in Fig. 
3. As can be seen, the ambient-pressure antiferromag- 
netic corundum- type phase, having a single kind of poly- 
hedral arrangement, is the most stable structure up to 
a pressure of about 34 GPa. Starting from this pres- 
sure the distorted corundum or Rh 2 03(II)-type phase in 
the LS ferromagnetic configuration is the thermodynami- 
cally stable phase. Several LS antiferromagnetic configu- 
rations are present a few mRy above the FM ground state 
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in the LS pressure range (see inset of Fig. 3), support- 
ing the conclusion that the system is above its magnetic 
ordering transition at the room temperature used in the 
MS experiments 22 . In our calculations, opv-type Fe203, 
often considered to be another possible candidate struc- 
ture for the high pressure phase, is never the most stable 
one in this pressure range. 
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FIG. 3: (Color online) Relative enthalpy curves of Fe203 
for corundum type, Rfi203(II) type and orthorhombic 
perovskite-type (opv) structures in different magnetic config- 
urations: ferromagnetic (FM), antiferromagnetic (AF), with 
HS or LS, and non-magnetic. Inset shows different LS AF 
configurations for Rh20s(II) type Fe203 with respect to the 
LS FM Rh 2 3 (II) type Fe 2 3 . 

All the magnetic configurations, in all investigated 
structures, show a spin collapse from ~ 3.8 /i#/Fe atom 
to ~ 1 ii b /Fe atom at some critical pressure. In 
hematite, the crossover from the HS to LS occurs at 
around 20, 25 and 36 GPa for the three different an- 
tiferromagnetic structures consistent with the elemental 
unit cell in Fig. 1(a), while for the ferromagnetic one 
it occurs at ~ 5 GPa. This value is ~ 30 GPa for 
the RI12 03(H) type Fe2 03 in the most stable magnetic 
(antiferromagnetic) configuration. Also opv-type Fe2 03 
displays a sharp spin transition as a function of pres- 
sure. In this structure there are two types of polyhedral 
cages (with 6-fold and 8-fold coordination) surrounding 
the cations and while Fe ions in 8-fold coordinated poly- 
hedra loose their spin gradually, Fe ions in the octahedral 
cages display a sudden drop from HS to LS at ~ 21 GPa. 
In all cases a volume collapse of about ~ 7 - 10 % at 
the spin transition is found, consistent with a first order 
transition and in agreement with experimental evidence. 
Hence, regardless of the considered crystal structures, the 
Fe atoms sitting in octahedral cages display a HS to LS 
transition under pressure. 

From the experimental point of view it is still contro- 
versial whether the spin transition occurs in the high- 
pressure or low-pressure structural phase. We believe 



that this spread in experimental results reflects the fact 
that spin transition is a general phenomenon occurring in 
all relevant structures and with a simple common origin 
as shown below. 



III. LOCAL STRUCTURE 

In order to clarify the interplay between spin and struc- 
tural transitions, we want to identify a few structural 
indicators that could be used to characterize the HS to 
LS transition and, as in previous studies, 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 we 
will concentrate on the Fe coordination sphere. Pressure 
and/or temperature induced (anisotropic) strain causes 
deformation of the octahedra, which can be quantified 
by parameters such as the Fe-0 bond-length distortion, 
SdF e -o — Y^h=i where D is the average octa- 

hedral Fe-0 bond-length and Adi is the deviation of the 
i-th Fe-0 bond- length from the average, and the octahe- 
dral angle variance, 34 given by AO 2 = Yl]=i(0j~ 90) 2 /ll, 
where 0j is the j-th O-Fe-0 angle whose ideal value is 
90°. 

In all structures, the octahedra become more regular at 
the transition as can be inferred from the sudden jump of 
both parameters, from higher to lower values, reported in 
Table I. However, Sdpe-o and AO 2 display very different 
values for different structures and hence these parameters 
alone are not suitable as a quantitative fingerprint of the 
occurrence of the spin transition. 

Other octahedral structural parameters can be defined 
such as the maximum (BL max ) and minimum (BL min ) 
Fe-0 bond- lengths, the octahedral volume (V oc t) an d 
a parameter measuring; the nearest-neighbor interaction, 
defined as 7^ = J2j=i e~ dj / r ° , where to = 1.3 A is 
of the order of the Oxygen ionic radius and dj's are the 
Fe-0 bond lengths. As can be seen in Table I, BL max ^ 
BL m i n and Y ct show jumps to lower values across the 
transition. Moreover they present very similar values in 
all structures, both just before but especially just after 
the transition when the spread of these parameters in 
different structures is of the order of a few percent. Sim- 
ilarly, the interaction parameter, 7^, also show a jump, 
this time from a lower to a larger value, again with a 
spread of a few percent. 

It is worth stressing that, in spite of the fact that the 
HS to LS transition occurs at different pressures in differ- 
ent structures, and the results accumulated contain both 
GGA and GGA -\-U calculations, for all the structural pa- 
rameters considered here, the spread in their calculated 
values, just before and especially just after the transition, 
is smaller than the jump they experience across the tran- 
sition, thus revealing the common mechanism operating 
in the different structures. All the above three structural 
parameters are therefore good indicators that could be 
used to monitor the approaching of a system to its criti- 
cal pressure. 

The octahedral volume is possibly the most funda- 
mental one and in Fig. 4 we report the evolution un- 



TABLE I: Calculated octahedral parameters at the GGA and GGA + U level: bond distortion index, 5dFe-o, octahedral angle 
variance, AO 2 , maximum bond length (BL macc ), minimum bond length (BL m ^ n ), interaction strength, and octahedral 

volume (Voct) at the HS — > LS transition. Three different antiferromagnetic configurations for the corundum- type phase are 
named as AF1, AF2 and AF3 and the most stable ground state magnetic (antiferro) configuration is considered for Rh 2 03(II) 
and opv structures. Overall averages, and their variation (in parenthesis), for all structural parameters are also indicated. Note 
that GGA + U calculations are done only for structures and magnetic configurations that are relevant (having relative stability 
at certain pressures) at the GGA level. And the values of U used here are our calculated values for the AF1 corundum type 
(3.3 eV) and FM Rh 2 3 (H) type (5.6 eV) Fe 2 3 . 
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FIG. 4: (Color online) Octahedral volumes for the HS and 
LS states. The data include corundum- type, Rh 2 03-type and 
optotype structures at the GGA level as well as for corundum- 
type and Rh 2 03-type structures at the GGA+U level. 



der pressure of the octahedral volumes of the different 
structural phases in the HS and LS states. The data 
include corundum- type, Rh203(II)-type and opv-type in 
the GGA scheme and, for corundum- type and RI12 03(H)- 
type, in the GGA+U scheme (£7=3.3 and 5.6 eV, respec- 
tively, as calculated according to Ref. 32). The tran- 
sition pressures 35 range from ~ GPa to ~ 70 GPa 
and the crystal packing is different for different struc- 



tures. In spite of all the above differences, the HS and LS 
octahedral volumes follow two well denned equations of 
state. The gap between these two equations of state cor- 
responds to the octahedral volume collapse at the transi- 
tion. This gap is larger at lower pressures and slowly de- 
creases with pressure. The finite width of these two equa- 
tion of states can be attributed to the fact that different 
structure display different inter-polyhedral surroundings, 
which are known to play some role in controlling the spin 
transition, 6 ' 7 as well as the fact that the results from dif- 
ferent DFT scheme (GGA and GGA+U) are collected 
together. Yet, it is remarkable that the finite width in 
the equation of states for the octahedral volume is only 
about 9% for the HS state and just ~ 3% in the LS one. 

A possible explanation for the larger uncertainty dis- 
played in the HS equation of state lays in the fact that 
in this state, owing to larger Fe-0 bonds and angular 
distortions (see table I), the influence of the surround- 
ings on the octahedra is more significant thus inducing 
larger spread in the octahedral volume. On the contrary, 
in the denser LS phase the octahedra become more regu- 
lar (table I) and the increased Fe-0 bond strength, 7^, 
reduces the relative influence of the surroundings, thus 
contributing to a smaller spread in the equation of state. 

Based on the above two regions in Fig. 4 it might 
also be possible to estimate the amount of volume col- 
lapse or the spin transition pressure at OK for similar 
kind of materials, provided any one of the quantities is 
known. Temperature effects on the transition pressure 
are however dramatic and are presently beyond our un- 
derstanding. 
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IV. DISCUSSION 

Given a reliable and consistent structural parameter 
to characterize the spin transition, let us now turn to an 
analysis of the main experimental facts involved in its in- 
terplay with the structural one. The co-existence of both 
low- and high-pressure structural phases over a finite 
width of pressure/temperature has already been pointed 
out by a number of studies 22 ' 24 ' 25 . Some controversy still 
exists 22 ' 23 on the nature of the intervening phases. The 
resistivity profile in [22] clearly identifies three pressure 
regions: a magnetic insulating phase at lower pressure, a 
nonmagnetic metallic phase at very high pressure and a 
mixture of them in between. In particular the nonmag- 
netic nature of the metallic phase in the intermediate re- 
gion was assigned on the basis of the assignment of a dou- 
blet in the Mossbauer spectra as arising from nonmag- 
netic Fe atoms. We notice however that this assignment 
is not completely compelling since also a para-magnetic 
phase, above its ordering transit ion- temperature, would 
contribute a doublet in MS. This interpretation would 
be consistent with the small energy spread between the 
different ferromagnetic configurations found in our cal- 
culations and reported in the inset of Fig. 3. Indeed, 
the x-ray K/3 emission spectra in [23] showed evidence of 
(para-)magnetism, although very weak. This study did 
not address the conductivity of the various phases. 

The picture arising from our calculations agrees with 
all the experimental evidence and helps clarifying the sit- 
uation. As can be seen from Fig. 3, the system has three 
broad regions in terms of magnetic states: a HS region 
at low pressure (R-T), a LS, but still (para-) magnetic, 
phase at intermediate pressure (R-II) and a nonmagnetic 
phase at very high pressure (R-III). While the HS phase is 
insulating, both corundum- type and Rli2 03(II)-type LS 
phases as well as the nonmagnetic phases are metallic. 
The HS insulating phase is followed by the LS metal- 
lic phase in either of the structural phases 37 , HP1 or 
HP 2, and we believe it is not constrained by any sharp 
boundary. In this crossover region, for a finite width on 
the pressure scale, the HS ^ LS transformation occurs 
in any particular octahedron depending on the volume 
that the particular octahedron acquires in that condi- 
tion. Our explanation for the spin transition goes as fol- 
lows: for a particular octahedron, the Fe-0 bond lengths 
and the octahedral volume decreases gradually with in- 
creasing pressure; at a certain point in the pressure scale, 
when the minimum Fe-0 bond length (see Table I) in 
that particular octahedron reaches a value for which the 
octahedral crystal field dominates over the exchange en- 
ergy, that particular octahedron acquires a new spin state 
switching from its 3 t2 g 2 e g state to b t2 g °e g , resulting in 
the collapse of the octahedral volume. Similar octahedra 
with equivalent surroundings will experience this transi- 
tion at the same pressure. While the transition pressure 
will vary for structures with different crystal packing. 

Our observation that HS ^ LS transformation depends 
on a local property, the octahedral volume, is in agree- 



ment with and explains the observed sluggishness of the 
structural transition 24 ' 25 . Inhomogeneity in the experi- 
mental pressure and temperature fields, the presence of 
defects and impurities in the experimental samples lo- 
cally affect the octahedral volume around Fe ions and 
therefore upon increasing pressure the HS to LS transi- 
tion will occur at different pressures in different octahe- 
dra depending on when the critical octahedral volume we 
have identified is reached. The local volume-collapse as- 
sociated with spin transition will contribute to enhance 
the inhomogeneity of the pressure-field in the sample and 
the transition will proceed gradually in a nucleation-and- 
growth way. 

In the same pressure range also the structural tran- 
sition occurs. Minor crystallographic changes, brought 
about by the alteration of one of the coordinating oxy- 
gen atoms, lead to distorted corundum type or Rli203(II) 
type from the corundum type structure. An important is- 
sue is whether spin and structural transitions are strongly 
inter-dependent or not. Our results indicate that the two 
type of transitions are not strongly dependent, rather the 
HS to LS transition is a general phenomenon occurring 
in all structural phases. This conclusion is in agreement 
with experimental evidence. Badro et. al 23 has indeed 
shown that the spin transition is not necessary for the 
structural transition to occur, and vice versa. Hence, 
based on our results and experimental findings, one can 
qualitatively conclude that the spin and structural tran- 
sitions are independent phenomena, in the case of Fe203, 
at least. 

The HS -> LS and HP1 -> HP2 transition are com- 
pleted by the upper bound of the intermediate region 
(R-II), which then progresses slowly to the nonmagnetic 
HP2 phase. Our calculated boundary for LS HP2 — > non- 
magnetic HP2 is ~ 120 GPa at zero temperature, which is 
rather far from the room temperature experimental value 
of ~ 72 GPa [22]. Even assuming an, expected, strong 
temperature effect the discrepancy remains large. Still, 
accepting the inadequacy of DFT in predicting the nu- 
merical value of the transition pressures, we stress that 
DFT at the GGA level predicts correctly the sequence 
and details of the transitions. Taken together, our DFT 
results and the two main experimental works 22 ' 23 provide 
a clearer picture of the pressure evolution of the system 
across three regions: a HS magnetic phase (R-I), a LS, 
but still (para-) magnetic, phase at intermediate pressure 
(R-II) and a nonmagnetic phase at very high pressure (R- 

iii). 

To further elucidate the evolution of the system across 
the spin transition and display the common mechanism 
that drives it in different crystal structures we present 
in Fig. 5 spin-density difference plots for HP2 (upper 
panels) and opv (lower panels). Before the spin transi- 
tion (see 1(a) and 2(a)) all Fe atoms are characterized 
by spherical spin-density difference contour plots, corre- 
sponding to the fact that the level occupation is 3 t2 g 2 e g 
for HS. Just after the spin transition (see 1(b) and 2(b)) 
Fe atoms in octahedral coordination sites switch to the 
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LS state. Their spin-density difference is characterized 
by four-lobe, due to the local 5 t2 g °e g level occupation. 
Notice that in panel 2(b) only the spin-density difference 
of the octahedral sites has changed nature; Fe atoms in 
the 8-fold coordinated sites remain HS. Finally panel 1(c) 
shows the spin density difference plot of HP2 when the 
system looses its spin completely and becomes nonmag- 
netic. On each site the averaged spin density vanish. The 
increased proximity of Fe and O atoms with increasing 
pressure enhances further the splitting of the d orbitals 
along with the pressure induced broadening of the states 
and a simultaneous transfer of some charge from O-p to 
Fe-d. This eventually leads to a spin state character- 
ized by 6 t2 g °e g level occupation. The similarity in the 
spin-density difference obtained in corresponding sites in 
different structures clearly show the common mechanism 
at work. 



V. CONCLUSIONS 



In conclusion, we have shown that, in spite of the gen- 
erally believed intricate relations with different structural 
parameters, the phenomenon of spin transition can well 
be described by a rather simple local structural parame- 
ter, ie., the octahedral volume, at least in Fe203, if not 
in general. The HS and LS octahedral volumes, separated 
by a well defined volume gap, are governed by two dis- 
tinct equations of state. For different structural phases, 
density functional theory, both at the GGA and GGA+U 
levels, reproduces this spin transition at a critical value 
of the octahedral volume. Analysis of our calculations 
suggests that at intermediate pressures the system is in 
a LS paramagnetic phase, thus clarifying some experi- 
mental controversy. 




FIG. 5: (Color online) Spin-density difference plots for the 
HP2 and opv across the HS — > LS transition. Upper panels, 
1(a), 1(b) and 1(c), refers the spin-density difference to HP2 
in the R-I, R-II, and R-III regions, respectively. Lower pan- 
els, 2(a) and 2(b), display the spin-density difference for opv 
before and after the HS — > LS transition, respectively. For the 
HP2 structure the plot is for (010) plane at 0.25 interplanar 
spacing away from the origin. For opv (100) plane at half of 
the interplanar spacing from the origin is shown. Note the 
spherical nature of the spin-density difference for Fe inside 
the 8-fold coordinated cage after the transition. The planes 
are antiferromagnetic, as can be seen in the color differences. 
For convenience the scale is shown in the left and two ex- 
treme points correspond to maximum values in the positive 
and negative direction. 
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